Rapidly rotating boson molecules with long or short range repulsion: 

an exact diagonalization study 
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The Hamiltonian for a small number, N < 11, of bosons in a rapidly rotating harmonic trap, 
interacting via a short range (contact potential) or a long range (Coulomb) interaction, is studied 
via an exact diagonalization in the lowest Landau level. Our analysis shows that, for both low and 
high fractional fillings, the bosons localize and form rotating boson molecules (RBMs) consisting of 
concentric polygonal rings. Focusing on systems with the number of trapped atoms sufficiently large 
to form multi-ring bosonic molecules, we find that, as a function of the rotational frequency and 
regardless of the type of repulsive interaction, the ground-state angular momenta grow in specific 
steps that coincide with the number of localized bosons on each concentric ring. Comparison of 
the conditional probability distributions (CPDs) for both interactions suggests that the degree of 
crystalline correlations appears to depend more on the fractional filling v than on the range of the 
interaction. The RBMs behave as nonrigid rotors, i.e., the concentric rings rotate independently of 
each other. At filling fractions v < 1/2, we observe well developed crystallinity in the CPDs (two- 
point correlation functions). For larger filling fractions v > 1/2, observation of similar molecular 
patterns requires consideration of even higher-order correlation functions. 

PACS numbers: 05.30.Jp, 03.75.Hh 
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I. INTRODUCTION 

Rotating ultra-cold trapped Bose condensed systems 
are most commonly discussed in the context of forma- 
tion of vortex lattices, which are solutions to the Gross- 
Pitaevskii (GP) mean-field equation 0, II 1 El S, IS 0, 
H, H, [l(J. Such vortex lattices have indeed been found 
experimentally for systems containing a large number 
of bosons [13, [ID, 131. Nevertheless, several theoret- 
ical investigations 1J, [ll| [H, E3, OH of rapidly ro- 
tating trapped bosonic systems suggested formation of 
strongly correlated exotic states which differ drastically 
from the aforementioned vortex-lattice states. While ex- 
perimental realizations of such strongly correlated states 
have not been reported yet, there is already a signifi- 
cant effort associated with two-dimenional (2D) exact- 
diagonalization (EXD) studies of a small number of par- 
ticles (N) in the lowest Landau level, LLL; the LLL 
restriction corresponds to the regime of rapid rotation, 
where the rotational frequency of the trap Q equals the 
frequency wn of the confining potential. The large ma- 
jority 0, [li, E3, El of such EXD studies have at- 
tempted to establish a close connection between rapidly 
rotating bosonic gases and the physics of elecrtons under 
fractional-quantum-Hall-effect (FQHE) conditions em- 
ploying the bosonic version of "quantum-liquid" analytic 
wave functions, such as the Laughlin wave functions, 
composite- fermion, Moore- Read, and Pfaffian functions. 

Recently, the "quantum-liquid" picture for a small 
number of trapped electrons in the FQHE regime has 
been challenged in a series of extensive studies [H, HF 



24j of electrons in 2D quantum dots un- 
der high magnetic fields (B). Such studies (both EXD 
and variational) revealed that, at least for finite sys- 
tems, the underlying physical picture governing the be- 



havior of strongly-correlated electrons is not that of a 
"quantum liquid." Instead, the appropriate descrip- 
tion is in terms of a "quantum crystal," with the lo- 
calized electrons arranged in polygonal concentric rings 
[Ll[2^j23^,[2l,[2l,[m[2^. These "crystalline" states 
lack 
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23j the familiar rigidity of a classical extended 
crystal, and are better described [H H3, El HI SI H 

as rotating electron ( or Wigner) molecules (REMs or 
RWMs). 

Motivated by the discovery in the case of electrons 
of REMs at high B (and from the fact that Wi gner 
molecules form also at zero magnetic field [H, |29|. l30l, 
l3ll . [321 . HH) some theoretical studies have most re- 
cently shown that analogous molecular patterns of lo- 
calized bosons do form in the case of a small number 
of particles inside a static or rotating harmonic trap 
[3~il [35l [H, [37], Hl| . In analogy with the electron case, 
the bosonic molecular structures can be referred to (36| 
as rotating boson molecules (RBMs); a description of 
RBMs via a variational wave function [33 ] built from 
symmetry-breaking displaced Gaussian orbitals with sub- 
sequent restoration of the rotational symmetry was pre- 
sented in Refs. HUE El. 



In this paper, we use exact diagonalization in the low- 
est Landau level to investigate the formation and prop- 
erties of RBMs focusing on a larger number of parti- 
cles than previously studied, in particular for sizes where 
multiple-ring formation can be expected based on our 
knowledge of the case of 2D electrons in high B. We 
study a finite number of particles (N < 11) at both 
low (y < 1/2) and high [y > 1/2) filling fractions 
v = N(N - 1)/2L [HJ (where L = C/h is the quantum 
number associated with the total angular momentum C) 
and investigate both the cases of a long (Coulomb) and 
a short ((^-function) range repulsive interaction. 
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As in the case of electrons in 2D quantum dots, we 
probe the crystalline nature of the bosonic ground states 
by calculating the full anisotropic two-point correlation 
function P(r, ro) associated with the exact wavefunction 
^(i"i, r 2 , rjy). The quantity P(r, r ) is proportional to 
the probability of finding a boson at r given that there 
is another boson at the observation point ro, and it is 
often referred to as the conditional probability distribu- 
tion [13, [H, [H |M HI H| (CPD). A main finding of 
our study is that consideration solely of the CPDs is 
not sufficient for the boson case at high fractional fill- 
ings v > 1/2; in this case, one needs to calculate even 
higher-order correlation functions, e.g., the full iV-point 
correlation function P(r; Ti,T2, rjv-i) [see Eq. ([12"]) ]. 

The present investigations are also motivated by re- 
cent experimental developments, e.g., the emergence of 
trapped ultra-cold gas assemblies featuring bosons inter- 
acting via a long-range dipole-dipole interaction. We 
expect the results of this paper to be directly relevant 
to such systems with a two-body repulsion intermediate 
between the Coulomb and the delta potentials. Addi- 
tionally, we note the appearance of promising experimen- 
tal techniques for measuring higher-order correlations in 
ultra-cold gases employing an atomic Hanbury Brown- 
Twiss scheme [4l[ or shot-noise interferometry [42|, EH . 
Experimental realization of few-boson rotating systems 
can be anticipated in the near future as a result of in- 
creasing sophistication of experiments involving periodic 
optical lattices co-rotating with the gas, which are ca- 
pable of holding a few atoms in each site. A natural 
first step in the study of such systems is the analysis of 
the physical properties of a few particles confined in a 
rotating trap with open boundary conditions (i.e., con- 
servation of the total angular momentum L) , such as we 
do here. 

Our main results can be summarized as follows: Simi- 
lar to the well-established (see, e.g., Refs. [22, 23]) forma- 
tion of rotating electron molecules in quantum dots, we 
describe here the emergence of rotating boson molecules 
in rotating traps. The RBMs are also organized in con- 
centric polygonal rings that rotate independently of each 
other, and the polygonal rings correspond to classical 
equilibrium configurations and/or their low-energy iso- 
mers. Furthermore, the degree of crystallinity increases 
gradually with larger angular momenta L's (smaller fill- 
ing fractions i/'s), as was the trend 0, [z| HH for the 
REMs and as was observed also for v < 1/2 in another 
study for rotating bosons in the LLL with smaller N 
and single-ring structures. We finally note that the crys- 
talline character of the RBMs appears to depend only 
weakly on the range of the repelling interaction, for both 
the low (see also Ref. [13]) and high (unlike Ref. (38J) 
fractional fillings j44j. 

In studies of 2 D q uantum dots, CPDs were used some 
time ago in Refs. [2(| [U [H, EH . For probing the intrinsic 
molecular structure in the case of ultracold bosons in 
2D traps, however, they were introduced only recently 
by Romanovsky et al. 34]. The importance of using 
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FIG. 1: Ground-state angular momenta, L gs , for N = 6 
bosons in a rapidly rotating trap [described by the LLL Hamil- 
tonian in Eq. ©] , as a function of the rotational frequency f2 
expressed in units of u>±. The bosons interact via a Coulom- 
bic repulsion and the many-body Hilbert space is restricted 
to the LLL. The angular momentum associated with the first 
bosonic Laughlin state occurs at L — 30, i.e., at N(N — 1). 
The value of c = 0.2&jxA for this plot; the many-body wave 
function does not depend on this choice. Note the stepwise 
variation of the values of the ground-state angular momenta, 
indicating the presence of an intrinsic point-group symmetry 
associated with the polygonal-ring structure of RBMs (see 
text for details). 



CPDs as a probe can hardly be underestimated. Indeed, 
while EXD calculations for bosons in the LLL have been 
reported earlier 0, HE [H, Ezl El] , the analysis in these 
studies did not include calculations of the CPDs, and 
consequently formation of rotating boson molecules was 
not recognized. 

This paper is organized as follows. In section [II] we es- 
tablish our notation and detail the physical assumptions 
that underlie the construction of our ab initio Hamilto- 
nian. In section IIII| we present our results for N = 6, 9 
and 11 bosons and compare the emerging crystallinity 
for both Coulomb and ^-function interactions. We sum- 
marize our results in section IV. The appendix presents 
explicit formulas for calculations of two-body matrix el- 
ements between two permanents. These permanents are 
used as the basis wave functions that span the many- 
boson Hilbert space; for the case of fermions, of course, 
one uses the more familiar Slater determinants 12011. 



II. THEORETICAL METHODS 

A. Ab initio Hamiltonian in the LLL 

Requiring a very strong confinement of the harmonic 
trap along the axis of rotation (huj z » Hu>±), freezes 
out the many body dynamics in the z-dimension, and 
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FIG. 2: Ground-state angular momenta, L gs , for N = 6 
bosons in a rapidly rotating trap [described by the LLL 
Hamiltonian in Eq. (J5)], as a function of the rotational fre- 
quency Q expressed in units of ui±. The bosons interact 
via a delta repulsion and the many-body Hilbert space is re- 
stricted to the LLL. The angular momentum associated with 
the first bosonic Laughlin state occurs at L = 30, i.e., at 
N(N - 1). The value of g = 2tvHlj±A 2 /N for this plot; the 
many-body wave function does not depend on this choice. 
The values of the angular momenta terminate with the value 
L — 30 (the Laughlin value) at Q/u>± = 1. In contrast, in 
the Coulomb-interaction case (see Fig. [TJ, the values of the 
ground-state angular momenta do not terminate, but diverge 
as fi/u>j_ - * 1- Note the stepwise variation of the values of 
the ground-state angular momenta in both cases, indicating 
the presence of an intrinsic point-group symmetry associated 
with the polygonal-ring structure of RBMs (see text for de- 
tails). 



the wavefunction along this direction can be assumed to 
be permanently in the corresponding oscillator ground 
state. We are thus left with an effectively 2D system. 
For such a setup, the Hamiltonian for N atoms of mass 
m in a harmonic trap rotating at angular frequency 
fiz is given by: 
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Here L = — HL = J2iLi^' r « x P« ^ s tne total angular 
momentum; r = (x,y) and p = {p x ,p y ) represent 
the single-particle position and angular momentum in the 
x — y plane and uj±_ is the radial trap frequency. 
The Hamiltonian can be rewritten in the form, 

H = f; ( (p * ~ m ™ x Tif + gfri - 



i=l 



2m 

N 



■^Vili-Tj). 



(2) 



The kinetic part of this Hamiltonian is formally equiva- 
lent to that of the Hamiltonian in the symmetric gauge of 
an electron (of charge e and mass m e ) moving in two di- 
mensions under a constant perpendicular magnetic field 
B, if one makes the identification that the cyclotron fre- 
quency ui c — eB/(m e c) — ► 2VL. 

We are now ready to examine our main approximation 
which is based on two assumptions: (1) that the rota- 
tional frequency fl is close to that of the confining trap; 
in this case the external confinement can be neglected in 
a first approximation and thus the single-particle spectra 
are organized in infinitely-degenerate Landau levels that 
are separated by an energy gap of 2fuv±, and (2) that the 
interaction strength is so weak that the mixing of Landau 
levels can be ignored. Since we work at zero temperature 
it then follows that all particles are in the lowest Landau 
level. This physical regime is known in the literature as 
the weak interaction limit. In a recent paper [46] it has 
been verified by explicit calculation that this approxima- 
tion, i.e., the lowest Landau level approximation, is very 
good for v < 1/2 (low fractional fillings), where strong 
crystalline correlations are most developed, as reported 
in the present work. 

We note, however, that strong correlations beyond the 
Gross-Pitaevskii mean field persist in the LLL even for 
v > 1/2 (high fractional fillings), as was found via EXD 
calculations for smaller N in Ref. [38[. In particular, 
Ref. HH concluded that even in the most favorable case 
(i.e., for L = N) formation of vortices for small N is 
not a prevalent phenomenon, and their appearance "is 
restricted to the vicinity of some critical values of the 
rotational frequency f2." This conclusion concurs with 
the results of Ref. [36| which used RBM variational wave 
functions to study a broad range of rotational frequencies 
that do not lead to the restriction to the LLL. The present 
work brings new insights to these findings by showing 
that the crystalline correlations for v > 1/2 (high frac- 
tional fillings) are clearly seen in the iV-point correlation 
functions, even though an inspection of the CPDs alone 
may be inconclusive (see examples in Section III. A). 

Taking into account that the single-particle spectrum 
associated with the Hamiltonian ^ (Fock-Darwin spec- 
trum HHH) is S iven b y e ™ = h \i 2n + \ l \ + !) W -L - ifi], 
the restriction to the LLL requires n = (Fock-Darwin 
single-particle states with zero radial nodes) and reduces 
the Hamiltonian H above [Eq. @] to the simpler form: 
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for which only the interaction term is non-trivial, since 
the many-body energy eigenstates are eigenstates of the 
total angular momentum as well; L = Yli=i 

In this study we shall investigate and compare our re- 
sults for strongly correlated states in the LLL obtained 
for a system of bosons interacting via a long or short 
range force respectively represented by the coulomb in- 
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FIG. 3: (a) Single-particle densities [n(r); left column], (b) CPDs [P(r, tq)] in 3D plots (middle column), and (c) CPDs in 
contour plots (right column), portraying the strengthening of the crystalline RBM structure for N = 6 bosons interacting via 
a repulsive Coulomb interaction as the filling fraction v is reduced. The white dots in the CPD plots indicate the reference 
point ro. We note in particular the gradual enhancement of the peak at the center of the plots, and the growth of the radius 
of the outer ring; the latter reflects the nonrigid-rotor nature of the RBMs (in analogy with the findings of Ref. [23j] regarding 
the properties of rotating electron molecules). The cases of v = 1/4 and v — 1/8 exhibit a clear (1, 5) crystalline arrangement, 
while the case of v = 1/2 (first Laughlin state) is intermediate between a (1,5) and a (0,6) pattern (see text for details). 
Lengths in units of A. The vertical scales are in arbitrary units, which however do not change for the panels within the same 
column (a), (b), or (c). 



ter action, 

«W = y y ( 4 ) 

or the contact potential, 

v(r) = gS 2 (r). (5) 

As may be inferred from the form of J/lll the cou- 
pling constants c and g trivially change the result by a 
multiplicative factor and their importance occurs only in 
comparing our results to a particular experimental sys- 
tem. Mathematically expressed, our description is valid 
only when: 

gA~ 2 < 2Hluj_ < hw z (6) 
in the case of a contact potential, and 

T < 2hu ± < haj z , (7) 



for the coulomb interaction, with A = ^h/(mu>±) be- 
ing the characteristic length of the harmonic trap. In 
these inequalities we have, for the sake of completeness, 
also included the energy scale associated with the 'frozen' 
perpendicular direction to give a full picture of the hier- 
archy of energy scales implicit in our approximations. 



B. EXD wave function and Conditional Probability 
Distribution 

For the solution of the large scale, but sparse, ma- 
trix eigenvalue problem associated with the Hamilton- 
ina Hlll, we have used the ARPACK computer code 
p^. [5fl|. For a given L, the EXD many-body wave func- 
tion is a linear superposition of permanents made out of 
the single-particle LLL wave functions, 
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FIG. 4: A ground state with a well-developed (0, 6) arrangement for N = 6 bosons interacting via a Coulomb repulsion at 
angular momentum L gs = 36. In classical electrostatic calculations [52J, the (0,6) arrangement is the first metastable isomer 
to the most stable (1,5) one shown by the CPDs of Fig. [3] The individual panels are: (a) 3-D plot of the single-particle density 
n(r); (b) 3D plot of the CPD; and (c) Contour plot of the CPD. In the CPD [P(r, To)] plots, the white dots indicate the 
reference point ro. Lengths in units of A. In (b) and (c), the vertical scales are the same. 
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with complex coordinates z = x + iy. Namely, 

K 



*(zi,z 2 ,...,zjv) = J2 Cjq> - 



(8) 



(9) 



,7=1 



where $j is a normalized permanent (see the Appendix); 
the index J denotes any set of N single-particle states 
{(j>i 1 (z 1 ),4> l2 (z 2 ), ...,4>i N (z N )} (not necessarily distinct) 
with angular momenta {h,h, ■ In} such that 



A' 



h < h < - < In, and ^l k = L. (10) 



fc=i 



The dimension K of the Hilbert space is controlled by 
the maximum allowed single-particle angular momentum 



^max, such that l k < Zmax, 1 < k < N. By varying 
l mn , we have checked that this procedure produces well 
converged numerical results. 

We probe the intrinsic structure of the EXD eigen- 
states [Eq. ^j] for crystalline behavior using the condi- 
tional probability distribution P(r,ro) defined (see, e.g., 
Ref. 3d]) by the following expression: 



N N 



P(r, r ) = <*| £ E 5 ^ - r W r * - r ^ 



(11) 



Unlike the case of electrons in QDs |22|], we found that, 
for rotating bosons, the CPDs are not sufficient to deci- 
pher the intrinsic molecular configuration in the case of 
high fractional fillings v > 1/2. In such a case, one needs 
to calculate higher correlation functions. In this paper, 
we use the iV-point correlation function defined as the 
modulus square of the full many-body EXD wave func- 
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FIG. 5: Contour plots of the CPD (a) and JV-point correlation 
function (b) for N = 6 bosons with L gs = 12 interacting via a 
Coulomb repulsion. The white squares indicate the positions 
of the fixed particles. The black square in (b) indicates the 
position of the 6th particle according to the classical (0, 6) 
molecular configuration. Note that the CPD in (a) fails to 
reveal the (0, 6) pattern, which, however, is clearly seen in 
the iV-point correlation function in (b). Lengths in units of 
A. The vertical scales are arbitrary. 
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FIG. 6: Contour plots of the CPD (a) and iV-point correlation 
function (b) for N = 6 bosons with L gs = 12 interacting via 
a 5-repulsion. The white squares indicate the positions of the 
fixed particles. The black square in (b) indicates the position 
of the 6th particle according to the classical (0, 6) molecular 
configuration. Note that the CPD in (a) fails to reveal the 
(0, 6) pattern, which, however, is clearly seen in the A-point 
correlation function in (b). Lengths in units of A. The vertical 
scales are arbitrary. 
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FIG. 7: Contour plots of the CPD (a) and iV-point correlation function (b) and (c) for N = 6 bosons with i gs = 15 interacting 
via a Coulomb repulsion. The white squares indicate the positions of the fixed particles. The black square in (b) and (c) indicates 
the position of the 6th particle according to the classical (1,5) molecular configuration. Note the different arrangements of the 
five fixed particles, i.e., (b) one fixed particle at the center and (c) no fixed particle at the center. Note also that the CPD in 
(a) fails to reveal the (1,5) pattern, which, however, is clearly seen in the iV-point correlation functions in both (b) and (c). 
Lengths in units of A. The vertical scales are arbitrary, but the same in (b) and (c). 
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FIG. 8: Contour plots of the CPD (a) and iV-point correlation function (b) and (c) for N = 6 bosons with L gs = 15 interacting 
via a 5-repulsion. The white squares indicate the positions of the fixed particles. The black square in (b) and (c) indicates 
the position of the 6th particle according to the classical (1, 5) molecular configuration. Note the different arrangements of the 
five fixed particles, i.e., (b) one fixed particle at the center and (c) no fixed particle at the center. Note also that the CPD in 
(a) fails to reveal the (1,5) pattern, which, however, is clearly seen in the iV-point correlation functions in both (b) and (c). 
Lengths in units of A. The vertical scales are arbitrary, but the same in (b) and (c). 



tion, i.e., 

P(r;ri,r 2 ,...,rjv-i) = |*(r; n, r 2 , r^-i)] 2 , (12) 

where one fixes the positions of iV — 1 particles and in- 
quiries about the (conditional) probability of finding the 
Nth particle at any position r. 

In the rest of this paper, we often find useful to contrast 
the CPDs and Appoint correlations to the single-particle 
densities n(r), which are circularly symmetric as a result 
of the total angular momentum being a good quantum 
number. The single-particle density is defined as 

N 

»(r) = <¥|X;*(ri-r)l*>- (13) 

i=l 



III. FORMATION OF ROTATING BOSON 
MOLECULES: NUMERICAL EXD RESULTS 

A. The N=6 case 

Since .Hlll is rotationally invariant, i.e., [H, L] = 0, its 
eigenstates must also be eigenstates of angular mo- 
mentum with eigenvalue hL. For a given rotational fre- 
quency f2, the eigenstate with lowest energy is the ground 
state; we denote the corresponding angular momentum 
as Lg S . 

We proceed to describe our results for N — 6 particles 
interacting via a Coulomb repulsion by referring to Fig. 
[1] where we plot L gs against the angular frequency Q. A 
main result from all our calculations is that L gs increases 
in characteristic (larger than unity) steps that take only 
a few integer values, i.e., for N — 6 the variations of L gs 



7 



COULOMB INTERACTION N=9 



DELTA INTERACTION N=9 



121— 

m 

107. 
100- 
93- 

R« 
79- 



72- 
64- 



56- 
51 — 



44. 
37 — 



„32- 



8—2- 



24 — 



0.65 0.70 NY 0.S 



.85 0.90 0.95 1.00 

FIG. 9: Ground-state angular momenta, L gs , for N = 9 
bosons in a rapidly rotating trap [described by the LLL Hamil- 
tonian in Eq. @] , as a function of the rotational frequency Q 
expressed in units of ui±. The bosons interact via a Coulom- 
bic repulsion and the many-body Hilbert space is restricted 
to the LLL. The angular momentum associated with the first 
bosonic Laughlin state occurs at L — 72, i.e., at N(N — 1). 
The value of c = 0.2huij_A for this plot; the many-body wave 
functions do not depend on this choice. Unlike the case of 
a 5-interaction (see Fig. I10[) . the values of the ground-state 
angular momenta do not terminate at L — 72 (first Laughlin 
state), but they diverge as fl/cu± — > 1. Note the stepwise 
variation of the values of the ground-state angular momenta, 
indicating the presence of an intrinsic point-group symmetry 
associated with the polygonal-ring structure of RMBs (see 
text for details). 



are in ste ps o f 5 or 6. In keeping with previous work on 
electrons [H, Hi HH, SI, H3, Lp at high B and very re- 
cently on bosons in rotating traps [34 l35l. l3ol l37l l38f . we 
explain these magic- angular-momenta patterns (i.e., for 
N = 6, L gs = Lq + 5k or L gs = Lq + 6k, with Lq = 0) as 
manifestation of an intrinsic point-group symmetry as- 
sociated with the many-body wave function. This point- 
group symmetry emerges from the formation of rotating 
boson molecules (RBMs), i.e., from the localization of 
the bosons at the vertices of concentric regular polygonal 
rings; it dictates that the angular momentum of a purely 
rotational state can only take values L gs = Lq + kirii, 
where rii is the number of localized particles on the ith 
polygonal ring [51]. Thus for N = 6 bosons, the series 
Lg S = 5fe is associated with an (1,5) polygonal ring struc- 
ture, while the series L gs = 6k relates to an (0, 6) arrange- 
ment of particles. It is interesting to note that in clas- 
sical calculations [HJ for N = 6 particles in a harmonic 
2D trap, the (1, 5) arrangement is found to be the global 
energy minimum, while the (0, 6) structure is the lowest 
metastable isomer. This fact is apparently reflected in 
the smaller weight of the L gs = 6k series compared to 
the Lgs = 5k series, and the gradual disappearance of 
the former with increasing L. 

Magic values dominate also the ground state angular 
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FIG. 10: Ground-state angular momenta, L gs , for N = 9 
bosons in a rapidly rotating trap [described by the LLL Hamil- 
tonian in Eq. ff3|> ] . as a function of the rotational frequency 
f2 expressed in units of u>±. The bosons interact via a delta 
repulsion and the many-body Hilbert space is restricted to 
the LLL. The angular momentum associated with the first 
bosonic Laughlin state occurs at L — 72, i.e., at N(N — 1). 
The value of g = 2irhujj_A' 2 /N for this plot; the many-body 
wave functions do not depend on this choice. Unlike the case 
of a Coulomb interaction (see Fig.O, the values of the ground- 
state angular momenta do terminate at L = 72 (first Laugh- 
lin state) when Q/lu± = 1. Note the stepwise variation of 
the values of the ground-state angular momenta, indicating 
the presence of an intrinsic point-group symmetry associated 
with the polygonal-ring structure of RMBs (see text for de- 
tails). 



momenta of neutral bosons (delta repulsion) in rotating 
traps, as shown for N = 6 bosons in Fig. [51 Although the 
corresponding f2-ranges along the horizontal axis may be 
different compared to the Coulomb case, the appearance 
of only the two series 5k and 6k is remarkable - pointing 
to the formation of RBMs with similar (1, 5) and (0, 6) 
structures in the case of a delta interaction as well (see 
also Refs. [13, HI]). An importance difference, however, is 
that for the delta interaction both series end at f2/ ui± = 1 
with the value L = N(N — 1) = 30 (the bosonic Laughlin 
value at v = 1/2), while for the Coulomb interaction this 
L value is reached for Q / loj_ < 1 - allowing for an infinite 
set of magic angular momenta [larger than ./V(7V — 1)] to 
develop as Q/u>± — > 1. Later, we will return back to this 
important difference between long range and short range 
interactions. 

Beyond the analysis of the ground-state spectra as a 
function of f2, the intrinsic crystalline point-group struc- 
ture can be revealed by an inspection of the CPDs [and 
to a much lesser extent by an inspection of single-particle 
densities] . Because the EXD many-body wave function is 
an eigenstate of the total angular momentum, the single- 
particle densities are circularly symmetric and can only 
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COULOMB INTERACTION N=9 (2,7) 1^=72+3x7 
(a) DENSITY (b) CPD1 (c) CPD2 




FIG. 11: Single-particle density and CPDs for the ground state of N = 9 bosons interacting via a Coulomb repulsion and having 
a total angular momentum of L gs = 93. The RBM crystalline structure corresponds to a (2, 7) polygonal-ring arrangement 
(which is the most stable classical equilibrium structure). Top row: contour plots. Bottom row: 3D plots. By column we 
describe the individual panels in both rows as follows: (a) Single-particle density, (b) CPD1 with reference point (indicated by 
a white dot) on the inner ring, revealing the localization of a second antipodal boson on the same inner ring. The outer ring 
exhibits a circular uniform symmetry that reflects the independent rotation of the two rings, (c) CPD2 with reference point on 
the outer ring revealing the localization of 6 bosons relative to the reference boson. We note that now the inner ring exhibits 
a uniform circular symmetry, reflecting again the independent rotation of the two rings. Lengths in units of A. The vertical 
scales are arbitrary, but the same in (b) and (c). 



reveal the presence of concentric rings through oscilla- 
tions in the radial direction. The localization of bosons 
within the same ring can only be revealed via the az- 
imuthal variations of the anisotropic CPD [Eq. (fTTj) ]. 
One of our findings is that for a given N the crystalline 
features in the CPDs develop slowly as L increases (or v 
decreases) . 

For v < 1/2, we find that the crystalline features are 
well developed for all sizes studied by us. In Fig. [3] and 
Fig. [U we present some concrete examples of CPDs from 
EXD calculations associated with the ground-states of 
A" = 6 bosons in a rotating trap interacting via a repul- 
sive Coulomb potential. In Fig. [31 we present the CPDs 
for L gs — 30 (bosonic Laughlin for v = 1/2), 60, and 
120; these angular momenta are associated with ground 
states at specific f2-ranges [see Fig. Q] . All three of these 
angular momenta are divisible by both 5 and 6. How- 
ever, only the L gs = 30 CPD (Fig. [3] top row) has a 
structure that is intermediate between the (1, 5) and the 
(0, 6) polygonal-ring arrangements. The two other CPDs, 
associated with the higher L gs — 60 and L gs = 120 ex- 
hibit clearly only the (1,5), illustrating our statement 
above that the quantum-mechanical CPDs conform to 
the structure of the most stable arrangement (i.e., the 
(1,5) for AT = 6) of classical point-like charges as the 
fractional filling decreases. On the other hand, in Fig. 
21 we plot the CPD for L gs = 36, which is divisible only 



by 6; we see that the corresponding single-particle den- 
sity and CPD are associated with a (0, 6) polygonal-ring 
arrangement. 

However, for v > 1/2, the azimuthal variations may 
not be visible in the CPDs, in spite of the characteris- 
tic step-like ground-state spectra [see Fig. [T] for A^ = 6 
bosons and Fig. [5] for N — 9 bosons]. This paradox is re- 
solved when one considers higher-order correlations, and 
in particular Af-point correlations [see Eq. (fl~2"]) ]. In Fig. 
[5] and Fig. [51 we plot the Af-point correlation functions 
for A" = 6 bosons and L gs = 12 for both the Coulomb in- 
teraction and (5-repulsion, respectively, and we compare 
them against the corresponding CPDs. The value of 12 
is divisible by 6, and one expects this state to be associ- 
ated with a (0, 6) molecular configuration. It is apparent 
that the CPDs fail to portray such sixfold azimuthal pat- 
tern. The (0, 6) pattern, however, is clear in the Af-point 
correlations, where we fixed five points (white dots) at 
the positions = ro exp(ifc-7r/3), k = 1,2,..., 5 (tq be- 
ing the radius of the maximum ring in the single-particle 
density) and we looked for the probability of finding the 
sixth boson at any other position r. The figures show that 
the maximum probability happens for r$ = ro exp(2ni) 
(black dot) which completes the (0, 6) regular polygon. 

A second example is given in Fig. [7] and Fig. [5] for 
the case of A" = 6 bosons and L gs = 15, and for both the 
Coulomb interaction and the (5-repulsion, respectively. In 
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COULOMB INTERACTION N=9 (1,8) 1^=72+2x8 
(a) DENSITY (b) CPD1 (c) CPD CONTOUR P 




FIG. 12: A ground state with a (1,8) RBM arrangement for N = 9 bosons interacting via a Coulomb repuslion and having 
angular momentum L gs = 88. In the classical calculation [S^], this arrangement is isomeric to the most stable (2,7) one 
exhibited by the CPDs in Fig. 1111 The individual panels are as follows: (a) 3D plot of the single-particle density, (b) 3D plot 
of the CPD. (c) Contour plot of the CPD. In the CPD plots, the white dots indicate the reference point ro. Lengths in units 
of A. The vertical scales are arbitrary, but the same in (b) and (c). 



COULOMB INTERACTION N=ll (3,8) Ig s =110+8 
(a) DENSITY n (b) CPD1 p (c) CPD2 p 

y 0.18 y 0.04 y 0.04 




FIG. 13: Single-particle density and CPDs for the ground state of N = 11 bosons interacting via a Coulomb repulsion and having 
a total angular momentum of L gs = 118. The RBM crystalline structure corresponds to a (3,8) polygonal-ring arrangement 
(which is the most stable classical equilibrium structure [52J). Top row: contour plots. Bottom row: 3D plots. By column we 
describe the individual panels in both rows as follows: (a) Single-particle density, (b) CPD with reference point (indicated by 
a white dot) on the inner ring, revealing the localization of two bosons on the same inner ring relative to the reference boson. 
The outer ring exhibits a circular uniform symmetry that reflects the independent rotation of the two rings, (c) CPD with 
reference point on the outer ring revealing the localization of 7 bosons relative to the reference boson. We note that now the 
inner ring exhibits a uniform circular symmetry, reflecting again the independent rotation of the two rings (see text for details). 
Lengths in units of A. The vertical scales are arbitrary, but the same in (b) and (c). 



this second case, one expects a (1,5) molecular pattern, 
which however is not visible in the CPDs. To reveal this 
(1,5) pattern, one needs to plot the TV-point correlations 
(middle and right panels). One now has two choices for 
choosing the positions of the first five particles (white 
dots), i.e., one choice places one white dot at the center 
and the other choice places all five white dots on the ver- 
tices of a regular pentagon. For both choices, as shown by 
the contour lines in the figures, the position of maximum 
probability for the sixth boson coincides with the point 



that completes the (1, 5) configuration [see the black dots 
at Tq = ro exp(27ri) in the middle panels and at the center 
in the right panels] . 



Note that for both cases, the differences in the CPDs 
and Appoint correlation functions between the Coulomb 
and the ^-repulsion are rather minimal. We will return 
to this result again in Section III.D. 
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COULOMB INTERACTION N=9 ^=72+7 DELTA INTERACTION N=9 1^=72+7 

(a) DENSITY (b) CPD (a) DENSITY (d) CPD 




FIG. 14: Comparison of single-particle densities and CPDs for three degenerate states (right half of figure) of N = 9 bosons 
interacting via a repulsive (^-interaction and having L gs = 79 against the lowest three (and non-degenerate) eigenstates (left 
half of figure) of the Coulomb repulsion for the same angular momentum. The numbers on the left hand side of the rows 
indicate the order of the eigenstates in the case of the Coulomb interaction, with 1 (top row) indicating the ground state, 
while 2 and 3 indicate excited (rovibrational 28]) states. The top-row plots reveal a (2, 7) crystalline (2, 7) structure for both 
types (long-range and short-range) interactions. Lengths in units of A. The vertical scales are arbitrary, but the same for the 
single-particle densities [(a) and (c)] and for the CPDs [(b) and (d)]. 



B. The TV = 9 case with Coulomb repulsion 

The results for N = 6 bosons portray many of the 
important features for all the iV's that we have exam- 
ined, keeping in mind that for N = 6 the corresponding 
polygonal-ring structures involve only a single nontriv- 
ial ring, i.e., they are of the form (1,5) or (0,6). From 
the case of confined electrons [22|, [§3j], it is known that 
the first nontrivial two-ring structure appears for N = 9. 
Our results presented here confirm that this is also the 
case for N = 9 bosons. 

In Fig. [9l we display the LLL ground-state angular mo- 
menta obtained from EXD calculations for N = 9 bosons 
as a funtion of fl/coo and for the case of a Coulombic re- 
pulsion. The remarkable feature again is the variation 
of L gs in well defined steps, with the steps taking only 
three specific values 7, 8, and 9 (with the step value of 
9 appearing only in two instances). The coresponding 
magic-angular-momenta values correspond to the series 
L gs = 2ki + 7k 2 (Sr), L gs = 8k («S 8 ), and L gs = 9k («S 9 ), 
associated with the polygonal-ring configurations (2,7), 
(1,8), and (0,9), respectively. Note that, as O/wo — > 1, 
only the S7 series survives, a behavior that is again re- 
lated to the fact that the (2, 7) structures is the most sta- 
ble isomer according to classical elecrostatic calculations 
[12], with the (1, 8) structure being the lowest metastable 
isomer. 

In Fig. 111! we present the single-particle density and 



CPD for N = 9 bosons and L gs — 93 (with ki = 1 and 
fca = 13). Naturally, the single-particle density is circu- 
larly symmetric, but it clearly protrays two concentric 
rings along the radial direction. Furthermore, the(2, 7) 
intrinsic azimuthal crystalline pattern is revealed in the 
two displayed CPDs, one with the observation point be- 
ing located on the inner ring (middle panels) and the 
other on the outer ring (right panels). It is remark- 
able that positioning the reference point ro on one ring 
reveals only the intrinsic crystalline correlations of the 
same ring, while the other ring appears uniform, and 
vice versa. This behavior is similar to that of multi- 
ring rotating electron molecules [12, [Hj] , and it suggests 
that the rings of the RBMs rotate independently from 
each other. It naturally leads to a physical picture of 
the RBMs as highly nonrigid rotors in analogy with the 
analysis [22, 23] of the rotating electron molecules in high 
B (see in particular sections V and VI in Ref. [23[). 



We close this subsection by showing a case associated 
the second isomer. Fig. [12] displays the single-particle 
density and CPDs for N = 9 bosons and L gs = 88 (be- 
longing to the 6>8 series with k = 11). Both of these quan- 
tities confirm that the ground state has a (1,8) intrinsic 
crystalline pattern, in agreement with the analysis above 
regarding the ground-state angular-momentum spectrum 
in Fig. M 
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(a) DENSITY 




DELTA INTERACTION N=ll (3,8) Lg^llO+8 
(b) CPD1 P 




(c) CPD2 




FIG. 15: Single-particle density and CPDs for one of the degenerate ground state of TV = 11 bosons interacting via a delta 
repulsion and having a total angular momentum of L gs = 118. The RBM crystalline structure corresponds to a (3, 8) polygonal- 
ring arrangement (which is the most stable classical equilibrium structure Top row: contour plots. Bottom row: 3D 
plots. By column we describe the individual panels in both rows as follows: (a) Single-particle density, (b) CPD with reference 
point (indicated by a white dot) on the inner ring, revealing the localization of two bosons on the same inner ring relative to 
the reference boson. The outer ring exhibits a circular uniform symmetry that reflects the independent rotation of the two 
rings, (c) CPD with reference point on the outer ring revealing the localization of 7 bosons relative to the reference boson. We 
note that now the inner ring exhibits a uniform circular symmetry, reflecting again the independent rotation of the two rings 
(see text for details). Lengths in units of A. The vertical scales are arbitrary, but the same in (b) and (c). 



C. The N=ll case with Coulomb repulsion 

Figure [T2] provides another example of the independent 
rotation of the rings that form the RBMs. In particular, 
Fig. 02] displays the single-particle density and the two 
CPDS (for the inner and outer ring) for L gs = 118 and 
N = 11 bosons interacting via a Coulomb repulsive po- 
tential. The intrinsic crystalline structure has a (3, 8) 
pattern, in agreement with the decomposition of the to- 
tal angular momentum as 118 = 3 x 2+8 x 14 and the fact 
that the (3, 8) isomer is the most stable one according to 
classical electrostatic calculations 1531 ■ 



D. Short-range vs. long-range interactions 

In the present investigation, we also also studied the 
influence of the range of the interaction on the properties 
of the many-body wavefunctions associated with a finite 
number of rapidly rotating bosons. In particular, using 
both a Coulombic and a delta interaction, we have calcu- 
lated and compared single-particle densities and CPDs 
for several sizes and fractional fillings v. Our result is 
that both types of interactions yield similar CPDs, i.e., 
that the degree of crystalline correlations does not de- 
pend strongly on the range of the interaction, as long as 
one refers to the same size N and fractional filling v. 

For example, the well developed crystallinity for v < 



1/2 familiar from the Coulombic case (see earlier sec- 
tions) maintains also for the case of a ^-interaction, as 
Fig. QH (N = 9) and Fig. [15] (N = 11) illustrate. In 
particular, Fig. [14] compares for N = 9 and L = 93 
the single-particle densities and CPDs associated with 
the three lowest states in the case of the Coulomb in- 
teraction (left half of the figure) against the correspond- 
ing quantities for three states out of the set of the now 
multiply degenerate [l5[ ground states in the case of the 
(^-interaction. Of course, this multiple degeneracy is un- 
physical, and in experimental applications one should use 
a short-range potential with a finite interaction length 
that is intermediate between the Coulomb and the 5- 
potential, like the dipole-dipole interaction [5J|. For our 
purposes here, it is sufficient to view the 5-interaction 
as a representative one of all short-range potentials, and 
thus to conclude that RBMs with similar patterns and 
properties do form, and that the degree of crystallinity 
is strengthened as v decreases, for both long-range and 
short-range interactions. We stress that the similarities 
involve both the crystalline arrangements [i.e., the (2, 7) 
polygonal-ring pattern] and the property that the rings 
rotate independently of each other, as can be seen by 
a direct comparison of the left and right halves of Fig. 
[TH We note that certain similarities (regarding the for- 
mation of RBMs in the LLL) between the long-range 
(Coulomb) and short-range (Gaussian interaction) cases 
have also been reported }37| most recently in the context 
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of EXD/LLL studies of one-ring RBMs (with A < 7 and 
v < 1/2). For similarities in the case v > 1/2, see Section 
III.A. 

Finally Fig. [T5] displays the CPD associated with one 
of the degenerate ground states for L = 118 and N = 11 
bosons interacting via a ^-potential. In keeping with our 
conclusions in the previous paragraph, we point out the 
clear formation of an RBM with a (3, 8) crystalline pat- 
tern; furthermore, this RBM also exhibits the property 
that the two rings rotate independently of each other. 



IV. CONCLUSIONS 

We presented results for the exact diagonalization of a 
small number A < 1 1 of rapidly rotating bosons interact- 
ing via a repulsive (^-function or Coulomb interaction. We 
focused on the weakly interacting limit where the assump- 
tion of confinement to the lowest Landau level is valid 
and studied in particular particle numbers which in the 
classical limit lead to polygonal crystalline configurations 
with at least two concentric polygonal rings. For all frac- 
tional fillings, we observed formation of rotating boson 
molecules, i.e., the appearance of crystalline correlations 
in close similarity to those found for confined two dimen- 
sional electron systems in high B [H [H, HE [H [H, Hi] , 
as well as for trapped strongly repelling bosons in the con- 
text of a recently developed trial wave function [H, H|| . 

Our main result is that for small filling fractions, i.e., 
v < 1/2, the molecular character is well developed and 
explicit in the CPDs; for larger fractions v > 1/2, the 
molecular character is clearly revealed through inspec- 
tion of higher-order correlation functions. The bosons 
organize into concentric polygonal rings which rotate in- 
dependently of each other, as was earlier found (221 . [23| 
for the case of rotating electron molecules in high B. 
The molecular structure is reflected in the ground-state 
spectra, i.e., as a function of the rotational frequency 
of the trap, the ground-state angular momenta vary in 
characteristic steps that coincide with the number of lo- 
calized bosons on each polygonal ring (5o| . Comparison 
of results for (5-function and Coulomb interactions also 
indicates that the emergence of the RBM is more depen- 
dent on the fractional filling v than on the range of the 
inter-particle interaction. 



APPENDIX A: MATHEMATICAL FORMALISM 
AND METHODS 

1. formalism 



Let 



ni!n 2 ! 



A! 



■^Ml^2u Pal (l)u Pa2 (2)...u PaN (N), 

( A1 ) 

represent one of the properly normalized permanents in 
Eq. © that form the basis of the many-body Hilbcrt 
space for the bosonic problem defined by the Hamilto- 
nian in Eq. ([3]). Such a permanent is constructed with 
M(< A) orthonormal single particle orbitals u ai {i) (LLL 
states in our case). The orbitals are indexed by a single 
integer on and (i) stands for all coordinates of particle 
i. While i — 1,2, . . . , N, the indices on are not necessar- 
ily distinct. The summation ^2 P runs over all distinct 
permutations of A particles among the M orbitals with 
occupations (multiplicities) {ni, n2, . ■ . , 

It is straightforward to show that the permanent de- 
fined in Eq. (lAlj) is normalized to unity. Indeed its norm 
is given by the following expression: 



•* A 1 AT 



nx\n 2 \ ■ ■ ■ n M \ 
A! 



i=l,N 

x E E U K (1)«0 01 C 1 )--- u k N (N)u QaN (A), 
P Q 



which upon evaluating the integrals yields 



(A2) 



ni\n 2 \ ■ ■ .njvf! 
A! 



ni!n 2 ! • . -n M ] - 



A! 



p Q 

E« 



ni!ri2! . . . njvj! A! 

A! n\\u2 \ ■ ■ ■ n jy! 



1 



where in the last line the summation equals the number of 
distinct ways of arranging A particles among M orbitals 
with the jth orbital having multiplicity nj. 



Interaction Matrix Elements 



Acknowledgments 

This research is supported by the U.S. D.O.E. (Grant 
No. FG05-86ER45234) and the NSF (Grant No. DMR- 
0205328). Calculations were performed at the Center for 
Computational Materials Science at Georgia Tech and 
at the National Energy Research Scientific Computing 
Center (NERSC). 



We make use of i.e., Eq. (jAljl above, which repre- 
sents a properly symmetrized and normalized wavefunc- 
tion to calculate the matrix elements for the many-body 
interaction term which we denote here as 

Vi =XS'- (A3) 

i<j 

At the outset we note that Vij = v{vi — r/) is symmetric 
with respect to the exchange of particles i and j and due 
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to the symmetry of bosonic wavefunctions, this property 
is carried over to the expression: 

($i(n,r 2) . . • , t n ), wy#B(ri,r 2 , . . . , r N )) (A4) 



In other words Eq. (| A4[) doesn't change value under ex- 
change of all the coordinates of particles i and j with any 
other pair. So 



{V l ) AB =J2(n,v^ B ) 



N(N-1) 



(3> A ,v 12 <S> B ) 



and the matrix element for the many-body Hamiltonian 
due to the interaction is reduced to 

(Vi)ab = N{N ~ 1 \ ®% v 12 $ B ) (A6) 
From our calculation of the norm in section fA II above. 



it should be obvious that matrix element defined by Eq. 
(|A6p can be non-zero if <&a and $ b are the same or differ 
in the occupations of two orbitals from each one; the case 
when &a and <& b differ in the occupations of one orbital 
from each permanent does not occur in the LLL. For all 
other cases the matrix elements in Eq. (|A6|) are zero. 
Because of the normalization factors introduced by the 
occupation numbers rij , the evaluation of Eq. (|A6|) breaks 
down into five cases, one for the diagonal elements and 
four for the off-diagonal elements. We shall take them in 
turn. 



3. Diagonal Elements 

Substituting the expression for the normalized wave- 
function Eq. (|Aip into (|A6|) above, we obtain 



(Vi) 



AA 



N(N — 1) 



/ n *. 

J i=l,N 



n\\ . ..n M l 
N) 



J2u* Pai (l)u* Pa2 (2)v 12 u Qai (l)u Qci2 (2) u Pa _(j)u Qaj (j) 

P,Q j=3,N 



N(N-l) f ni \n 2 \...n M \ 
dridV2 Nl 



^^ qi (1) U ^ 2 (2) Ui2MQqi (1) UQq2 (2) 5 PaitQo 
P,Q 



i=3,N 



r 



(A7) 



In Eq. (|A7|) only two terms in the sum J^q can be 
nonzero, i.e., the terms satisfying the two conditions be- 
low 

Pai = Qa 1 and Pa 2 = Qa 2 i Or 
Poli — Qa 2 and Pa 2 — Qan 

and for both of the above cases 



Thus Eq. (|A7|) is rewritten as 



, N. 



(A8) 



nil . . . n a \np\ . 



,n M i 



(Vi)aa- 2(7V _ 2) , 
x^((a/3\v 12 \aP)+(ap\v 12 \0a)) ]T (1), 



a/3 



P{N-2} 



(A9) 



where we have simplified the notation of the indices, i.e., 
we set P ai — a and P a2 = /3, and we used the following 
shorthand definition for the matrix elements of the two- 
body interaction 



case to be considered. When a = (3, one has 

(N — 2)! 



E (D 



P{N-2} 

When a (3, one has 



E « 

P{N-2} 



(N-2)\ 



ml... (n a - 1)1(710 - 1)! . . .n M ! 



(All) 



(A12) 



With the above results, Eq. (|A9[) finally yields for the 
diagonal elements of the interaction part of the Hamilto- 
nian 



(Vi) 



AA 



n a (n a - 1) 



M 



~y^ i {aa\v\ 2 \aa) 



i a np 



M 



((aP\v 12 \aP) + (a(3\v 12 \(3a)) 



a7 tf3=l 



(A13) 



(a/3\v 12 \jA) = J dr 1 £fr 2 ti*(l)«J(2)«i2U 7 (l)«A(2). 

(A10) 

In Eq. (|A9|) , we have explicitly noted the fact that the 
sum refers now to N — 2 particles. There are two 



4. Off-Diagonal Elements 

Here, we extend the formalism developed in the previ- 
ous sections to be able to treat matrix elements between 



14 



two different basis permanents. Let and $^1/3-1 
denote permaments having the occupation numbers 

{ni,n 2 , . . .,na,np,n*y,nA, . . ■ 7 n M } (A14) 

and 

{ni,n 2j . . . ,^-1,71/3-1,7^+1,^+1, ■ ■ ■ ,n M }, (A15) 

respectively. Then we need to calculate the matrix ele- 
ment 



^a— l'fl-i can ^ e m t er P re t e d as the permanent created 
from when two particles are taken from the single 
particle states labelled by a and and placed into the 
states 7 and A. Four cases arise; they are 



7 + A 

7 ^ A 

7 = A 

7 = A 



a ^ 

a = 

a ? 

a = 



, $T +1,A+1 x, W 



(A16) 



J 



• Case I: 7 / A; a / 

We begin from Eq. (|A6[) and make adjustments for the difference in the occupation numbers to obtain: 

/ <5 *7+l,A+l y$ 7A\ _ 
V V a-l,/3-l i K I^a/3 J ~ 



N(N - 1) _ /m! ...{n a - 1)1(71/3 - !)!(% + l)!(riA + 1)! • • • n M \ _ /n i! . . . n a !n^!n 7 !n A ! ■ ■ ■ n M ! 

JV! 
(JV-2)! 



2 V AT! 

x(( 7 A|?j 12 |a/3) + <7A|w 12 |/3a) + (A 7 |v 12 |a^) + (A 7 |v 12 |/?a)) 



ni! . . . (n a - l)!(n/3 - l)!n 7 !?i A ! • . .n M ! 



or finally 



,^*7 + l,A+l y-*7A\ 



^n a np(n 7 + l)(n A + l)(<7^l«ia|a/3> + (7A|«ia|/3a}) 

With similar considerations we obtain the results for the other cases. 
Case II: 7 ^ A; a = 0. 



($*' 



* 7 +l,A+l y^jA 



N(N - 1) /m! ...(«„- 2)!(n 7 + l)!(n A + 1)! . . • n M ! /ni! . . . n Q !n 7 !n A ! . . • n M ! 



AH 



AH 



x ((7A|7Ji 2 |aa) + (A7|wi 2 |aa)) 



(N — 2)! 



ni! . . . (n a — 2)!n 7 !riA! • ■ ■ njst ! 



(A17) 



(A18) 



(A19) 



(A20) 



or finally • Case IV: 7 = A; a = /3. 

(C 7 - + 2 1,A+1 , = ( $ a T -2j ^2) = 

\]n a {n a - l)(n 7 + l)(n A + l)(7A|wi 2 \aa) -^\J n a (n a - l)(n 7 + l)(n 7 + 2)(77|t>i 2 |aa) 

(A21) (A23) 



• Case III: 7 = A; a ^ 0. 



^JnaUfsirij + l)(n 7 + 2)(77|?Ji 2 |a!^) 



5. Conditional probabilty distributions 

If the many-body wave function is denoted as 
V?(rx, r2, . . . , rjy), the conditional probability of finding 
(A22) a particle at position r given that there is another one at 
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position Tq is given by 

P(r,r ) = ^^ f J] dr,|vl/(r,r ;r3,...,r,v)| 2 . 

i=3,iV 

(A24) 

We reduce the calculation of the conditional probabil- 
ity (within a proportionality constant) to evaluating the 
expectation value of the following two-body operator 

T = Y / S(r-r l )S(r -r 3 ). (A25) 
If r ^ r , the expression under the sum is not symmet- 



ric with respect to interchanging rv and r, [unlike the case 
of the two-body potential i?(r,j — Vj)]. We rewrite T as: 

T = ^(JCr-rO^ro-r^+XJ^r-rO^ro-r,-) 

= £X<J(r - r,)<5(r - r,-) + *(r - r^fo - r<)). 

i<j 

(A26) 

which has a symmetric expression under the summation. 
Then the calculation of the CPD proceeds exactly as the 
calculation for the interaction part of the Hamiltonian. 
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